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ABSTRACT 


The  error  due  to  image  quantization  in  stereoscopically 
evaluated  range  is  analyzed,  using  a  simple  camera  model  and 
assuming  that  the  image  planes  of  the  two  cameras  are  coplanar . 
The  "best”  position  and  the  uncertainty  space  for  a  point  in 
space  are  calculated  on  the  basis  of  the  quantized  image  co¬ 
ordinates  of  that  point  in  the  two  cameras. 


ii 


1. 


Introduction 


In  this  technical  note  we  develop  the  mathematical  formulas  for  de¬ 
termination  of  range  values  and  their  extreme  errors  on  the  basis  of  quan¬ 

tized  stereoscopic  pictures  of  a  scene. 

Two  cameras,  with  lenses  L„  and  L  ,  are  located  at  a  distance  D  from 

Z  r 

each  other  to  determine  the  components  of  a  distance  r  from  the 

midpoint  between  L„  and  L  and  an  object  point  B,  as  shown  in  Fig.  1.  The 
Z  r 

focal  length  of  each  camera  is  f.  The  picture  plane  is  of  size  w  x  w  and  is 

quantized  into  n  x  n  cells,  each  of  size  d.  We  wish  to  calculate  the  error 

introduced  by  the  digitization  process  in  determining  the  x,  y,  and  z  of  r. 

Let  us  denote  the  coordinates  of  the  image  point  on  the  left  and  right 

pictures  by  (xfl,  z  )  and  (x  ,  z  ),  respectively.  For  convenience,  we  shall 
Z  Z  r  r 

use  mirrored  image  planes,  instead  of  the  actual  ones.  To  simplify  the 

analysis,  let  us  first  assume  that  z  =  z  =0,  i.e.,  the  projecting  rays 

Z  r 

(of  lengths  r  and  r  )  lie  in  the  plane  defined  by  the  parallel  optical 

Xj  T 

axes  of  the  two  cameras  (see  Fig.  1).  We  shall  later  remove  this  assump¬ 
tion.  Furthermore,  we  shall  also  assume  that  the  optical  axes  are  parallel 
to  the  y  axis;  the  effects  of  panning  and  tilting  of  the  cameras  may  be 
added  later  by  using  standard  transformation  rules. 

2.  Quantization  Errors  and  Uncertainty  Regions 

Thetop  view  of  the  stereopsis  plane,  in  which  both  x  >  0  and  x  >  0, 

Z  r 

is  shown  in  Fig.  2.  Due  to  the  quantization  process,  we  obtain  picture  co¬ 
ordinates  x^*  and  xr*>  and  thus  arrive  at  Point  B*  whose  coordinates  are  x* 
and  y*  instead  of  the  real  point,  Point  B,  whose  coordinates  are  x  and  y. 

By  inspection, 

y*  _  JL 

x*+^D  x^* 

and 

y*  f 

x*-£D  x  * 
r 


1 


(1) 


(2) 


For  given  quantized  values  of  x  *  and  x  *,  the  quantization  errors 

Xj  r 

6  and  6  (see  Fig.  2)  introduce  uncertainty  into  the  computed  values  of 

Xj  r 

x  =  x*  -  Ax  and  y  =  y*  -  Ay.  A  quadrilateral  uncertainty  region  is  ob¬ 
tained  by  varying  6„  and  6  according  to  -d/2  <  6/,<  d/2  and  -d/2  <  6  <  d/2, 

£  r  Xj  r 

as  shown  in  Fig.  3.  The  extreme  x  and  y  values  of  the  quadrilateral  un¬ 
certainty  region  determine  a  rectangle  defined  by  x  £  x  ^  x  and 

min  max 

y  £  y  £  y  ,  This  rectangle  constitutes  a  worst-case  uncertainty 
min  max  - 

region.  It  is  marked  in  Fig.  3  by  a  dotted  rectangle. 

The  values  of  6  and  6  in  the  above  four  cases  of  extreme  quantiza- 

Xj  r 

tion  errors  are  given  in  Table  I.  It  can  be  shown  that  the  values  of 

-6-6  and  -6+6  affect  the  expressions  for  x,  y,  Ax,  and  Ay.  The  re- 
£  r  £  r 

suiting  expressions  are  given  in  Table  I  for  each  of  the  four  cases.  Note 
that  in  Cases  2  and  3  (see,  for  example,  Fig.  3),  Ay  =  0  and  the  values 
of  Ax  are  equal  in  magnitude  but  opposite  in  sign. 

3.  Object  Zones 

So  far  we  have  assumed  that  x  *  >  0  and  x  *  >  0,  i.e.,  Point  B*  is 

£  r 

to  the  right  of  L  .  There  are,  however,  two  more  possibilities  to  consider: 
r 

x  *  >  0  and  x  *  <  0,  in  which  Point  B*  is  to  the  right  of  L.  and  to  the  left 

of  L  ,  and  x  *  <  0  and  x  *  <  0.  in  which  Point  B*  is  to  the  left  of  L„. 
r  £  r  £ 

Note  that  the  fourth  possibility,  x  *  <  0  and  x  *  >  0,  is  forbidden.  Note 

£  r 

also  that  by  excluding  x  *  =  0  and/or  x  *  =  0,  we  have  tacitly  assumed 

Xj  I* 

that  picture-cell  boundaries  exist  at  =  0  and  at  x^  =  0.  Generally, 
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Table  I.  Calculated  Expressions  for  x,  y,  Ax,  and  Ay  in  Four  Cases  of  Extreme  Quantization  Errors 


therefore,  there  are  three  possible  object  zones ,  as  shown  in  Fig-.  4: 

Zone  R  (x  *>0  and  x  *>0),  Zone  M  (x  *X)  and  x  *<0) ,  and  Zone  L  (x  *<0 
£  r  £  r  £ 

and  x  *<0).  Any  of  these  zones  is  characterized  by  the  algebraic  in- 
r 

equality  x^£  x^  (if  Point  B  is  at  infinity,  then  x^=  x^) .  Therefore, 
assuming  that  the  two  cameras  have  identical  quantization  cells,  then 
x  *  a  x  *. 

I  r 

Quadrilateral  and  worst-case  rectangular  uncertainty  regions  in  each 

of  the  above  zones  are  illustrated  in  Fig.  4.  Each  uncertainty  region  is 

defined  by  two  uncertainty  beams,  -d/2  <  6  „<  d/2  and  -d/2  <  6  <  d/2.  The 
-  £  r 

quadrilateral  (or  rectangular)  uncertainty  region  in  Zone  R  is  similar  to 

the  one  in  Fig.  3.  Its  extreme  boundaries  are  x  ,  =  x  ,  x  =  x  ,  y 

min  1  max  4  min 

=  yn ,  and  y  =  y  .  where  each  subscript  denotes  the  "case"  number,  as 
1  max  4  - 

defined  in  Table  I.  Keeping  the  positive  x  *  fixed,  we  "fan"  the  right 

Xj 

uncertainty  beam  counterclockwise ,  thus  causing  the  uncertainty  region  to 

move  to  the  left  and  into  Zone  M,  until  some  negative  x  *  is  reached. 

r 

Unlike  Zone  R,  the  extreme  x  boundaries  of  the  uncertainty  region  in  Zone 

M  are  x  =  x  and  x  =  x  ,  whereas  the  extreme  y  boundaries  remain  y 

min  2  max  3  min 

=  y  and  y  =  y  .  Next,  we  keep  the  negative  x  *  fixed  and  fan  the  left 

uncertainty  beam  counterclockwise,  again  causing  the  uncertainty  region 

to  move  to  the  left  until  reaching  the  position  in  Zone  L.  The  extreme 

x  values  have  changed  again,  i.e.,  x  =  x  and  x  =  x  .  and  the  extreme 

rain  4  raax  1 

y  values  are  again  y  =  y  and  y  =  y  . 

min  1  max  4 

4.  Quantization  Coordinate  Numbers 


So  far  we  have  used  the  quantization  coordinates  x*  and  x*  to  de¬ 
termine  the  values  of  x,  y,  and  their  quantization  errors.  To  be  more 


general,  we  may,  alternatively,  use  the  picture  indices  corresponding  to 
x*  and  x*.  Following  the  assumption  that  picture-cell  boundaries  exist 

jC  r 

at  x.=  0  and  x  =  0,  then 
£  r 


V  =  (V4)d 
v  =  <vJ)d  - 


n  .  £  n 
£  r 


(3) 


where  n  and  n 
i  r 


are  quantization  coordinate  numbers, 


which  are  nonzero 
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integers,  and  the  inequality  ns  n  follows  from  the  inequality  x  *  £  x  *. 

Hr  Hr 

The  quantized  image  shift  is  defined  as 


x  *  =  x  *  -x  *  =  (n  -n  ) d  (4) 

A  4  r  H  r 

Assuming  that  D,  f ,  and  d  are  constant,  we  note  in  Table  I  that  y  and  Ay 
in  every  case  and  Ax  in  Cases  2  and  3  are  each  a  function  of  the  image 
shift  only.  However,  x  in  every  case  and  Ax  in  Cases  1  and  4 

are  each  a  function  of  both  the  image  shift  and  the  average  quantized  x 
coordinate 


x  * 
P 


x  *+x  *  (n  +n  -l)d 
H  r  <L  r 


2  2 


(5) 


We  consider  first  y  and  Ay  vs.  the  quantized  image  shift,  n  -n  . 

i  r 

Based  on  the  expressions  given  in  Table  I,  plots  of  y  =  y,  and  y  =  y, 

rain  1  max  4 

vs.  n  -n  are  shown  in  Fig.  5.  As  expected,  both  y  and  y  increase 
i  r  min  max 


as  n  -n  decreases.  Note  that  y 


=  y 


max,n  -n  +2 
l  - 


Also  note  in 


l  "r  -  ’  ’  -----  —  ■'min.n  -n 

JL  T  Xj  t 

Fig.  5  that  as  n  -n  decreases,  the  worst-case  error  y  -y  increases. 

&  r  max  min 

More  information  about  the  error  in  y  may  be  obtained  by  separating  y 

max 

-y  into  its  components  Ay  =  Ay  and  Ay  =  -Ay  •  Based  on  Table  I, 
min  min  1  max  4 

plots  of  Ay  .  and  Ay  vs.  n  -n  are  shown  in  Fig.  6. 
mm  max  i  r 

The  errors  in  x  in  Zone  M  are  inversely  proportional  to  n^-n  >  whereas 

all  other  errors  (errors  in  x  in  Zones  R  and  L  and  errors  in  y  in  all 

2 

zones)  are  approximately  inversely  proportional  to  (n^-n^)  .  Furthermore, 
since  n  2  1  and  n  s  1,  then  n  -n  =  |n  |+|n  |  and,  therefore,  the  magnitudes 


a 


4 


of  the  errors  in  both  x  and  y  in  Zone  M  are  smaller  than  those  in  Zones  R 

and  L  for  the  same  x  *  and  x  *  magnitudes.  If  x  *  =  -x  *,  then  x  *  =  0, 

4  r  i  r  p 

i.e.,  Point  B*  is  along  the  axis  of  symmetry  between  the  two  cameras 
(cf.  Figs.  3  and  4).  For  this  special  case,  the  expressions  for  x  and  Ax 
in  Cases  2  and  3  (see  Table  I)  are  reduced  to 

D 


x  =  Ax  =  -x  =  -Ax 
max  max  min  min 


4n 


(6) 


i 
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Referring  to  Fig.  4,  the  area  of  the  uncertainty  region  in  Zone  M  is 
minimized  for  a  given  projecting  ray  if  the  two  projecting  rays  are  per¬ 
pendicular.  The  locus  of  such  minimum  regions  is  a  semicircle  centered 
at  x  =  0  whose  diameter  is  D.  Unfortunately,  this  locus  is  likely  to  be 
too  close  to  the  cameras  to  have  any  practical  value. 

5.  Vertical  Coordinates 


So  far  we  have  assumed  that  the  object  point  lies  on  the  x-y  plane. 

We  now  lift  this  restriction  and  let  Point  B  be  in  a  three-dimensional 
space,  as  shown  in  Fig.  7.  The  ideally  parallel  optical  axes  of  the  two 
cameras  point  along  the  y  direction,  as  before.  The  left  and  right  pro¬ 
jecting  rays  pierce  their  corresponding  image  planes  at  two  points  whose 

picture  coordinates  are  (x  ;  z  )  and  (x  ;  z  )  respectively.  The  above 

Z  Z  r  r 

analysis  for  determination  of  the  x  and  y  coordinates  assuming  z  =  0 

is  valid  with  no  change  for  the  general  case  in  which  z  /  0  because  x 

and  y  depend  only  on  x „  and  x  ,  not  on  z  „  or  z  . 

Z  r  Z  r 

If  there  were  no  difference  between  the  tilt  angles  of  the  two  cameras, 
the  two  projecting  rays  would  meet  at  Point  B  (see  Fig.  7).  In  this  case, 
the  intersection  of  the  plane  defined  by  the  two  rays  and  the  vertical 
plane  y  =  f  is  a  straight  line  parallel  to  the  x-y  plane  and  at  a  dis¬ 
tance  z  =  z  from  this  plane.  By  inspection  of  similar  triangles  in  Fig.  7, 
Z  r 

z/z  =  z  /z  =  f/y;  hence,  z  =  z  y/f,  and  the  "best  estimate"  of  z  is  z* 

Z  r  Z 

=  z  *y*/f.  Substituting  y*/f  =  D/(x  *-x  *)  from  (2),  and  replacing  z  * 

Z  Z  r  Z 

by  the  average  value 

z  *  =  £(z  *+z  *)  (7) 

p  Z  r 


just  in  case  z  *  /  z  *  due  to  error  in  camera  alignment,  we  find  that 
Z  r 


z*  -  Dz  */x  * 
P  A 


(8) 


Since  each  of  the  quantized  coordinates 
error  that  varies  between  -d/2  and  d/2, 


x  *,  x  *,  z  *,  and  z  *  may  have  an 
Z  r  z  r 

the  value  of  z  may  vary  between 


6 


z  and  z  ,  where 
min  max 

z  *-£d 

z  =  D  -P 
min  x  *+d 

A 


and 


z  *+|d 

z  =  D  — - 

max  x  *-d 

A 


(9) 


(10) 


6.  Three-Dimensional  Uncertainty  Space 

The  error  in  estimating  x,  y,  and  z  due  to  the  picture  quantization 

is  confined  to  a  hexahedral  uncertainty  space  formed  by  the  two  three- 

dimensional  uncertainty  beams  as  shown  in  Fig.  8  (note  that  x^*>D/2  in 

Fig.  8  by  mistake).  Each  of  the  uncertainty  beams  has  the  shape  of  a 

four-sided  pyramid.  The  side  views  of  these  beams  along  the  -x  axis  are 

shown  in  Fig,  8(a).  These  side  views  overlap  because  we  have  decided  to 

approximate  the  ranges  z  =  z*  +  4d  and  z  =  z*  +  ^d  by  the  average  range 

SL  i  ~  r  r  — 

£(z^*+zr*)+£d,  (Note  that  z^  or  z^  is  not  necessarily  within  this  average 
range.)  On  the  basis  of  this  approximation,  the  uncertainty  space  is  hexa 
hedral  with  four  vertical  sides,  as  shown  in  Fig.  8(b).  Without  this 
assumption  the  two  uncertainty  beams  may  not  meet  and  the  averaging 
process  will  become  unnecessarily  complex. 

Let  us  summarize  the  above  analysis  in  terms  of  the  quantization 
numbers 


n£x 

(X£*/ d) 

+ 

4 

(ID 

n  = 

(x  */ d) 

+ 

i 

(12) 

rx 

r 

n4z  = 

(z  */d) 

Hj 

+ 

i 

(13) 

and 

n  =  (z  */d)  +  J  .  '  (14) 

rz  r 

Using  the  abbreviations 
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and 


n  =  n  -n  =  x  */d 
A  Ax  rx  A 

n  =  §(n  +  n  )  =  (x  */d)  +  £ 

x  jlx  rx  p 


n  =  § (n  +  n  )  =  (z  */d)  + 
z  Az  rz  p 


(15) 

(16) 

(17) 


(4),  (5),  and  (7)  become 


and 


x  *  =  n  d 
A  A 


x  *  =  (n  -J)d 
P  x 


z  *  =  (n  -^)  d 
P  z 


(18) 

(19) 

(20) 


Hence,  (1),  (2),  and  (8)  are  expressed  as 


( 


v*  - 


x=M  I  D(n  -\)/n  ) 
x  A 


z* 

v  J  \ 


Df/(n  d) 

A 

D(n  -J)/n  , 

z  Ay 


(21) 


and  the  expressions  in  Table  I  in  conjunction  with  Fig.  4  are  reduced  to 
the  following: 


D(n  -£)/(n  +1)  if  n „  >0  and  n  >  0  (Zone  R) 
x  A  Ax  rx 


min 


=  {  D (n  -l)/n 


x 


if  n„  >  0  and  n  <0  (Zone  M) 
Ax  rx 


D(n  -£)/(n  -1)  if  n„  <  0  and  n  <0  (Zone  L) 
x  '  A  Ax  rx 


(22) 


x 


D(n  -|)/(n  -1)  if  n  >0  and  n  >  0  (Zone  R) 
x  A  Ax  rx 


max 


Dn  /n 
x  A 


if  n  >  0  and  n  <0  (Zone  M) 
Ax  rx 


D(n  -l)/(n  +1)  if  n  „  <  0  and  n  <  0  (Zone  L) 
x  A  Ax  rx 

Df 


min  d(n  +1) 
A 


(23) 


(24) 
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Df 


max  d(n  -1) 
A 


(25) 


min 


[  D(n  -l)/(n  +1) 
\  z  A 

LD(n  -l)/(n  -1) 
z  A 


if  n  >  | 


if  n  *  i 


(26) 


and 


max 


|'Dn  /(n  -1) 
z  A 


Dn  /(n  +1) 
z  A 


if  nM 


if  n  <  i 
z 


(27) 


The  relative  errors  in  x,  y,  and  z  are  found  from  (21)  through  (27) 
to  be  as  follows: 


<^> 

x 


X* 


1/n 


if  n  >  0  and  n  >0  (Zone  R) 
Xx  rx 


x 

min  min 


-1  =<  l/[2(n  -1)]  if  n  >  0  and  n  <  0  (Zone  M)  (28) 


i: 


l/n 


Xx 


rx 


if  n„  <  0  and  n  <  0  (Zone  L) 
Xx  rx 


(A-) 

X 


X* 


max  max 


-1/n 


-1  =<(  -1/ (2n  ) 

X 


l/n 


if  n  >  0  and  n  >  0  (Zone  R) 

Xx  rx 

if  n  >0  and  n  <0  (Zone  M)  (29) 
Xx  rx 

if  nA  <  0  and  n  <  0  (Zone  L) 

Xx  rx 


<^> 

y  y  .  n 

mm  mm  A 

(AZ)  =  JCL  -1  .  -  i 

y  y  n 

max  max  A 


(30) 


(31) 


(A5) 

z 


z* 


-1  = 


mm  mm 


(2n  +  n  -l)/[2(n  -l)nA]  if  n  >  \ 
z  A  z  A  z 

(~2n  +  n  +  1 )/[2 (n  -l)n  1  if  n  5  i 
z  A  z  A  z 


(32) 


(45) 

z 


z* 


-1  = 


max  max 


(“2n  -n  +l)/(2n  n  ) 
z  A  z  A 

(n  -n  -l)/(n  n  ) 
z  A  z  A 


if  n^ 
z 

if  n  <  | 
z 


(33) 


9 


7. 


Space  of  View 


The  effect  of  the  size  of  the  image  plane  will  now  be  considered, 
assuming  that  the  image  plane  is  a  w  X  w  square.  The  widest  space  viewed 
from  each  camera  is  a  pyramid  formed  by  the  straight  lines  connecting  each 
lens  with  the  corners  of  its  square  image  plane.  The  top  view  along  the 
-z  direction  and  the  side  view  along  the  -x  direction  of  the  two  pyramids 
are  shown  in  Fig.  9.  The  half-infinite  space  in  front  of  the  two  cameras 
is  divided  into  three  spaces:  (1)  the  space  not  seen  by  any  camera,  (2) 
the  space  seen  by  only  one  camera,  and  (3)  the  space  seen  by  both  cameras. 

The  top  and  side  views  of  the  third  space  are  shaded  in  Fig.  9,  and  each 
is  characterized  by  the  view  angle 

?  ='  2  tan  1  —  .  (34) 

2f 

Clearly,  only  this  tent-shaped  space  is  effective  for  stereo  applications. 
The  distance  between  the  origin  (midpoint  between  the  lenses)  and  the  tip 
of  the  effective  space  is 

r  =  f—  .  (35) 

m  w  v 

This  is  the  minimum  distance  that  can  be  determined  by  the  stereoscopic 
setup. 

8.  Numerical  Example 

Let  us  consider  a  numerical  example  in  which  D  =  6M  =  152  mm,  f  =  16  mm, 

w  =  10  mm,  and  the  image  plane  is  quantized  into  120  X  120  cells.  We  wish 

to  determine  the  effective  stereopsis  space,  the  analog  and  quantized  image 

3  4 

coordinates  associated  with  an  object  point  located  at  v  =  (10  mm,  10  mm, 

3  t 

2  •  10  mm)  ,  and  the  uncertainty  space  that  corresponds  to  the  image  quan¬ 
tized  data  (if  these  data  were  given). 

(1)  Effective  Stereopsis  Space 

Using  (34)  and  (35), 


10 


and 


6  =  w  tan  1  ~~  =  34.7° 
2*16 


152 

r  -  16  — —  =  243  ram. 
m  10 


(2)  Image  Coordinates 


The  quantization  cell  size  is 


d  = 


w 


10  1 

=  —  mm  =  0.0833  mm. 


120  120  12 


Modifying  (8) ,  z  = 


Df 

152*16 

y 

o 

r—i 

II 

D  X 

0.243 

152 

XA 

0.243 

z 

D 

152 

V 

we  get 

=  0.243  mm. 


10  -  1,6  mm. 


2  •  10  =  3.2  mm. 


x  =  h  (2x  +  x  )  =  i  (3.2  +  0.243)  =  1.7215  mm 
Z  p  A 


x  =  \  (2x  -  x  )  =  \  (3.2  -  0.243)  =  1.4785  mm 
r  P  A 


The  image  quantization  numbers  are 


V  rVd1*=  1 

1.7215/0.08331 

=  [20.7]  =  21 

n  =  |x  /d]  =  | 
rx  r 

1.4785/0.0833] 

=  ^7.71  =  18 

n  =  \k  /d]  =  | 
z  p 

[3.2/0.0833]  = 

[38.41  =  39 

Following  (11)  through  (20), 


x  *  =  (n  -  *)d  =  20.5  •  0.0833  =  1.71  mm 

i  JLx 

x  *  =  (n  -  A)d  =  17.5  •  0.0833  =  1.46  mm 
r  rx 


*  We  shall  use  the  notation 
M  =  round  up  integer  £x 
[ xj  s  round  down  integer  ^x 
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21  -  18  =  3 


V 


x  * 
P 

z  * 
p 


£  (21  +  18)  =  19.5 

3  •  0.0833  =  0.25  mm 

(19.5  -  0.5)  •  0.0833  =  1.585  mm 

(39  -  0.5)  •  0.0833  =  3.21  mm. 


Following  (21) , 


X* 

152 

*  19/3 

1  963  1 

V*  = 

y* 

=■ 

152 

*  16/(3*0.0833) 

= 

9,730 

i  z*  1 

152 

•  38 

.5/3  j 

1 1,950 

absolute 

errors 

are,  therefore, 

X* 

-  X 

1-  37 

Av  = 

y* 

-  y 

= 

-270 

(mm) 

z* 

"  Z! 

1  -50  i 

The  corresponding  relative  errors  are 


/  Ax/x 

j(x*-x)/x\ 

-37/1, 000 

Ay/y 

= 

(y*-y)/y 

- 

-270/10,000 

- 

Az/z 

(z*-z)/z 

-50/2,000 

(3)  Uncertainty  Space 


Following  (22)  through  (27), 


x  .  =  152  *  19/4  =  722.5  mm 

min 

x  -  152  •  19/2  =  1,445  mm 
max 

152-16 

y  =  « — — — ; —  =  7 , 300  mm 
min  0.0833-4 


y  =  152-16  „ „ 

max  — —  —  =  14,600  mm 
0.0833*2 


z  =  152  •  38/4  =  1,445  mm 
min 


z  =  152  *  39/2  =  2,962  mm. 
max 


(mm) 


-0.037 

-0.027 

-0.025 
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Thus,  the  best  estimate  and  the  uncertainty  space  of  v  =  (x,  y,  z)  are 

x*  =963  mm;  722.5  <  x  <  1,445  mm 
y*  =  9,730  mm;  7,300  <  y  <  14,600  mm 
z*  =  1,950  mm;  1,445  <  z  <  2,962  mm. 

The  resulting  relative  uncertainty  space  is  defined  fay  -0.33  <  &x/x  <  0.33, 

-0.33  <  Ay/y  <  0.33,  and  -0.34  <  Az/z  <  0.35.  These  results  agree  with 

(28)  through  (33)  in  which  n  =  3  and  n  =  39. 

A  z 
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